Molecular structures in charmonium spectrum: The XYZ puzzle 
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We study in the framework of a constituent quark model the possible contributions of molecular 
structures to the XYZ charmonium like states. We analyze simultaneously the cc structures and 
the possible molecular components in a formalism which allows us to treat channels below and above 
thresholds. The only molecular state found in the 1 ++ sector correspond to the X(3872). Molecular 
resonances also appear with other quantum numbers. So, the so called Y(3940) and the V(3915) are 
suggested to be J PC = ++ charmonium states. In the J PC = 1 sector we also found significant 
contributions of the molecular structures which can affect the phenomenology. 
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I. INTRODUCTION. 

In the last few years several charmonium like states 
were observed with similar masses near 3.9 GeV but with 
quite different properties and in very different production 
processes. Altogether these states were called the XYZ 
states. A complete list of these new states can be found 
in Ref. Among them we will only comment on the 
confirmed states. 

This new charmonium era started around 2003 when 
the Belle Collaboration discovered the lightest one, the 
X(3872), in the exclusive decay B ± -> K^+ir- J/tp 0. 
The mass of the state was measured to be 3872.0 ± 0.6 
MeV very close to the M D o +M D *o threshold. The width 
was found to be very small T < 2.3 MeV. The state 
was soon confirmed by CDF [|, DO and BaBar @. 
By combination of the recent results reported by the 
Belle @, BaBar Q and CDF Q Collaborations, the mass 
value is established at M x = 3871.55 ± 0.20 MeV. 

An striking feature of the A (3872), which cannot be 
explained by a simple cc structure, is the ratio Q 



Ri 



V(3872) 



-7r°J/> 



V(3872) ->■ n+ir-J/iP 



0.8 ±0.3. (1) 



The dipion mass spectrum in the ir + Tr~J/ip channel 
shows that the pions come from the p° resonance. On 
the other hand the 7r + 7r~7r° mass spectrum has a strong 
peak around 750 MeV suggesting that the process is dom- 
inated by an uj meson. Although this number should be 
corrected by the strong phase suppression of the ujJ/rp 
channel against the pJ/ip one, the ratio R\ ~ 1 is incom- 
patible with a traditional charmonium assumption and 
is telling us that some isospin mixing is needed and that 
this mixing requires the contribution of both neutral and 
charged DD* channels. Recently the Belle Collabora- 
tion [10l | measured the ratio 



Ri = 



r(V(3872) -> 7^(25)) 
T(A(3872) -> 7J/VO 



< 2.1(at 90% C.L.) (2) 



which complicates the interpretation of this state. 

Two years later, the so called at that time V(3940) 
was observed by the Belle Collaboration as a near- 



threshold enhancement in the coJ/ip invariant mass dis- 
tribution for the B — > KuiJ/ip decay [ll[. Belle re- 
ported a mass of M = 3943 ± 11 ± 13 MeV and a width 
r = 87±22±26 MeV. Belle observation seems to be con- 
firmed by BaBar [3, although the mass (M = 3914±4.1 
MeV) and the width (r = 29±10 MeV) were both smaller 
than Belle values. A later measurement of BaBar [l3| 
confirmed this last mass value (M = 3919 ± 3.8 ± 2.0 
MeV). Very recently a new charmonium-like state, the 
V(3915), has been reported by the Belle Collaboration 
in the 77 — > J/i/jui decay [lif . The measured mass is 
M = 3914 ± 3 ± 2 MeV and the width r = 23 ± 9 MeV. 
It has not yet been seen in the DD channel. Despite of 
the different mass and width of the first measurement of 
Belle, some authors [l|, [l5[ subsumed both states under 
the name A(3915), although the question if there is one 
or two different resonances is not definitely settled. 

Another charmonium like state, the V(3940), was ob- 
served in this region by Belle as a resonance in the double 
charmonium production e + e~ — > J/ipDD* in the mass 
spectrum recoil against the J/ip [16| . Later on Belle con- 
firmed the observation of A(3940) ->■ DD* decay [l7j . 
In addition Belle found a new charmonium like state 
V(4160) decaying into D*D* . Neither of them have been 
seen in the experimentally more accessible DD channel. 
As the decay of the Y73940) — > DD* was not observed in 
the B Y(3940)K pj], the A(3940) and the V(3940) 
should be different states. 

Two more states increase the experimental findings in 
this region. The Z(3930) was reported by Belle in the 
DD channel produced in 77 collisions with mass and 
width M = 3929 ± 6 MeV and T = 29 ± 10 MeV 0. 

Finally, it is worth to mention that one more resonance 
has been classified as 'well established' in Ref. [l| . It was 
found by the BaBar Collaboration [2(| in the reaction 
e+e" -> DD with a mass of 3943 ± 17 ± 12 MeV and a 
width of 52 ± 8 ± 7 MeV and confirmed by Belle [2l|. It 
was called G(3900) in Ref. 

Concerning the quantum numbers of the new states, 
measurement of the angular correlations between final 
state particles in the A(3872) — > Tr + ir~J/?p decay [22| 
together with small phase space available for the decay 
A(3872) ->■ D D tt° [H strongly favors the J PC = 1++ 



2 



quantum numbers for this state. However some recent 
results seems to favor negative parity for this meson [l3| . 

The situation is worse in the case of the AT (3940). 
It has not been seen in the DD channel which rules 
out the J p = ++ assignment. The dominant DD* 
mode suggest that the AT (3940) is the cc(2 3 Pi) state with 
J PC — 1 ++ but this quantum numbers seems to be as- 
signed to the X(3872). 

More consensus exists with the assignment of the 
Z(3930). The DD decay mode makes it impossible to be 
the 77 c (35) J PC = state. The two photon produc- 
tion can only produce DD states in ++ or 2 ++ and these 
two cases can be distinguished looking to the dN/dcosO 
distribution being 6 the angle between the incoming pho- 
ton and the D meson in the 77 center of mass system. 
This distribution is flat for ++ states and behaves like 
sin 4 8 for 2 ++ . The measurement of Belle strongly 
favors the 2 ++ case. 

The J PC assignment in the case that the Y"(3940) was 
a cc state is still unclear. A conventional cc charmonium 
interpretation is in principle disfavored since it is well 
above the threshold for open charm decays and then these 
decay modes would dominate over the <j)J/ip and ujj/rp 
decay rates. The X(3915) is preliminary assigned to be 







++ 



2++ \U 



Finally due to the entrance channel of the production 
reaction the G(3900) is clearly a J PC — 1~~ state. 

Obviously it is difficult to accommodate all these states 
in a qq scheme and all type of hypothesis about their 
structures (molecules, hybrids, tetraquarks) has been 
proposed in the literature (see [24[ for a review). 

The P = +, C = + sector is specially suited for the 
coexistence of cc states and molecular structures. The 
reason is the following: taken into account the negative 
intrinsic parity of the quark-antiquark pair, to get a pos- 
itive parity state one needs at least one unit of angular 
momentum. However four quarks can reach the same 
positive parity with zero angular momentum. Then, the 
energy increase due to the angular momentum excitation 
may compensate the two additional light quark masses 
making the cc and the ccqq structures almost degenerate. 

This mechanism has been suggested in 25j as a possi- 
ble explanation to the AT (3872) properties. In this refer- 
ence a coupled channel calculation of the 1 ++ sector has 
been performed including cc and DD* molecular config- 
urations. The AT (3872) appears as a dynamically gener- 
ated cc and DD* molecule by the coupling to a Xci (2-P) 
quark state. Although the cc mixture is less than 10%, it 
is important to bind the molecule. In addition tt + tt~ J/ip 
decay modes data from Belle and BaBar are reasonably 
explained. 

In this paper we will propose a theoretical explanation 
of the nature of some of these states using the same con- 
stituent quark model of Ref. [25[ . This has been success- 
fully used to describe hadronic spectroscopy and hadronic 
reactions [23429] ] and is based on the assumption that 
constituent quark mass is a consequence of the sponta- 
neous chiral symmetry breaking and has been recently 



applied to the study of the energy spectrum and d ecay 
properties of the J PC = 1 charmonium sector [30]. 
Taken into account the existence of several thresholds in 
this energy region we develop a formalism which treats 
simultaneously molecular states above and below the dif- 
ferent thresholds. This allows us to consider both dynam- 
ically generated molecular states by the coupling with cc 
structures and molecular components in the dressed cc 
states. 

The paper is organized as follow. In Section [II] we first 
discuss the basic ingredients of the constituent model, the 
coupled channel formalism and the coupling mechanism 
between the different channels. Results and comments 
are given 111 Section [nil Finally we summarize the main 
achievements of our calculation in Section [IVI 



II. THE MODEL 
A. The constituent quark model 

The constituent quark model used in this work has 
been extensively described elsewhere [281 ] and therefore 
we will only summarize here its most relevant aspects. 
The chiral symmetry of the original QCD Lagrangian 
appears spontaneously broken in nature and, as a conse- 
quence, light quarks acquire a dynamical mass. The sim- 
plest Lagrangian invariant under chiral rotations must 
therefore contain chiral fields, and can be expressed as 



C = 4>{ip- M(q 2 )U^)ip 



(3) 



where U lb 



e 7- 



-4> a i5 



is the Goldstone boson fields ma- 



trix and M{q 2 ) the dynamical (constituent) mass. This 
Lagrangian has been derived in Ref. [3lJ as the low- 
energy limit in the instanton liquid model. In this model 
the dynamical mass vanishes at large momenta and it 
is frozen at low momenta, for a value around 300 MeV. 
Similar results have also been obtained in lattice calcula- 
tions [32[ • To simulate this behavior we parametrize the 
dynamical mass as M(q 2 ) = m q F(q 2 ), where m q ~ 300 
MeV, and 



F{q 2 ) 



A 2 

v 



A 2 

x 



r 



(4) 



The cut-off A x fixes the chiral symmetry breaking scale. 

The Goldstone boson field matrix U 15 can be expanded 
in terms of boson fields, 



C/ 75 = 1 + - r 7 5 A°7T a - 7-^TT a TT 



(5) 



The first term of the expansion generates the constituent 
quark mass while the second gives rise to a one-boson ex- 
change interaction between quarks. The main contribu- 
tion of the third term comes from the two-pion exchange 
which has been simulated by means of a scalar exchange 
potential. 
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In the heavy quark sector chiral symmetry is explicitly 
broken and this type of interaction does not act. However 
it constrains the model parameters through the light me- 
son phenomenology and provides a natural way to incor- 
porate the pion exchange interaction in the open charm 
dynamics. 

Below the chiral symmetry breaking scale quarks still 
interact through gluon exchanges described by the La- 
grangian 



C gqq = iy/4na s 1p~1aG^\ c lj) , 



(6) 



where A c are the SU(3) color generators and G£f the 
gluon held. The other QCD nonperturbative effect corre- 
sponds to confinement, which prevents from having col- 
ored hadrons. Such a term can be physically interpreted 
in a picture in which the quark and the antiquark are 
linked by a one-dimensional color flux-tube. The sponta- 
neous creation of light-quark pairs may give rise at same 
scale to a breakup of the color flux-tube [33|. This can 
be translated into a screened potential [34| in such a way 
that the potential saturates at the same interquark dis- 
tance 

Vcon (n 3 ) = {-oc (1 - e~"« r » ) + A} (A* c < • X°j ) . (7) 

Explicit expressions for these interactions are given in 
Ref. El. 



B. The coupled channel calculation 

In this section we present the formalism for the cou- 
pling of molecular structures with the cc spectrum. We 
start defining the meson wave functions we will use all 
along the paper. To found the quark-antiquark bound 
states we solve the Schro ding er equation using the Gaus- 
sian Expansion Method f36|. In this method the radial 
wave functions solution of the Schrodinger equation are 
expanded in terms of basis functions 



(8) 



where a refers to the channel quantum numbers. The co- 
efficients b™ and the eigenenergy E are determined from 
the Rayleigh-Ritz variational principle 



E 



(9) 



where the operators T", n and N% n are diagonal and the 
only operator which mix the different channels is the po- 
tential V$n- 

A crucial problem of the variational methods is how 
to choose the radial functions 4>nii r ) m order to have a 
minimal, but enough, number of basis functions. Follow- 
ing (36j we employ gaussian trial functions whose ranges 
are in geometric progression. The geometric progression 



is useful in optimizing the ranges with a small number of 
free parameters. Moreover the distribution of the gaus- 
sian ranges in geometric progression is dense at small 
ranges, which is well suited for making the wave function 
correlate with short range potentials. The fast damp- 
ing of the gaussian tail is not a real problem since we can 
choose the maximal range much longer than the hadronic 
size. 

To model the cc system we assume that the hadronic 
state is 

|*> =^^0.) +Y,MP)\M B fi) (10) 

a 

where \ip a ) are cc eigenstates of the two body Hamil- 
tonian, are qq eigenstates describing the A and B 
mesons, \(Pa<Pb/3) is the two meson state with j3 quan- 
tum numbers coupled to total J PC quantum numbers 
and x/3 (P) is the relative wave function between the two 
mesons in the molecule. When we solve the four body 
problem we also use the gaussian expansion of the qq 
wave functions obtained from the solution of the two 
body problem. This procedure allows us to introduce in a 
variational way possible distortions of the two body wave 
function within the molecule. To derive the meson-meson 
interaction from the qq interaction we use the Resonating 
Group Method (RGM). 

The coupling between the two sectors requires the cre- 
ation of a light quark pair nn. Similar to the strong de- 
cay process this coupling should be in principle driven by 
the same interquark Hamiltonian which determines the 
spectrum. However Acklch et al. 13711 have shown that 
the quark pair creation 3 Po model [38|], gives similar re- 
sults to the microscopic calculation. The model assumes 
that the pair creation Hamiltonian is 



H = g (fx ip(x)ip(x) 



(11) 



which in the non-relativistic reduction is equivalent to 
the transition operator [39] 



T 



3 V2 7 ' J2 J d3 P d3 P' 6(3) (P + P') 



yi 



p-p 



bUpHdp 1 ) 



C=1,/=0,S=1,J=0 



(12) 



where /x (y = fl) are the quark (antiquark) quantum num- 
bers and 7' = 2 5 / 2 7r 1 / 2 7 with 7 = s 2 - is a dimensionless 
constant that gives the strength of the qq pair creation 
from the vacuum. From this operator we define the tran- 
sition potential hp a (P) within the 3 Pq model as 

/3|T|V a > =Php a (P)6W(P cm ) (13) 

where P is the relative momentum of the two meson 
state. 

Adding the coupling with charmonium states we end- 
up with the coupled-channel equations 
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c a M a + J2 I h a p(P) X p(P)P 2 dP = Ec a 
J2 J H p ,p(P',P) X0 (P)P 2 dP + J2h0, a (P')c a = Exp'(P') 



(14) 



r 



where M a are the masses of the bare cc mesons and Hp> p tained from the qq interaction. Solving the coupling with 
is the RGM Hamiltonian for the two meson states ob- the cc states we arrive to an Schrodinger-type equation 

I 



£ J (H f} , fi {P>.P) + V;ff{P\P))x P {P)P 2 dP = E X p l (P') 



(15) 



r 



where 



V${P',P;E) 



h > a (P')h aP (P) 



E — M n 



(16) 



Our aim is to find molecular states above and bellow 
thresholds in the same formalism. However, above the 
threshold we will find complex eigenenergies, where the 
imaginary part is related to the width of such states. 



J 



In order to find the poles of the T matrix we must be 
in the correct Riemann sheet, so we have to analyti- 
cally continue all the potentials for complex momenta. 
Once the analytical continuation is performed, the previ- 
ous coupled channel equations can be solved through the 
T{p,p ,E) matrix, solution of the Lippmann-Schwinger 
equation, 



Tf > 'f ) {P',P;E) = V$' fi {P',P ] E) + Y J J V^" (P', P"; E) — ^ T^(P", P; E) P" 2 dP 



(17) 



where V^' P {P',P;E) = V^'^(P', P) + V^//(P', P; E), 

V 'P(P',P) is the RGM potential and Vpf/ (P' , P; E) is 
the effective potential due to the coupling to intermediate 
cc states. 

In this way we study the influence of the cc states on 
the dynamics of the two meson states. This is a differ- 
ent point of view from the usually found in the literature 



r 



where the influence of two meson states (in general with- 
out meson-meson interaction) in the mass and width of 
cc states is studied 40]. Our approach allows to gener- 
ate new states through the meson-meson interaction due 
to the coupling with cc states and to the underlying qq 
interaction. 

The T matrix of Eq. (fTTj) can be factorized as [4l[ 



(18) 



T^(P>, P- E) = T^(P', PE) + Y, 4> P ' a ' (P'; E)A~} a (E)r (P; E) 

i 

with A -1 (E) = ((E — M )8 a ' a + Q aa (E)) * being ^ ma trix of the RGM potential excluding the coupling 

^ / to the cc pairs. 

the propagator of the mixed state and T{, (P', P; E) the „, c . iBn/n m u • , j ^ 

° ^ 1 he new functions qr (P; E) can be interpreted as the 
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dressed 3 Pq vertex by the RGM meson-meson interaction 
and are defined as 



(P;E) 



Tf{P,q-E)hp, a (q) 
q 2 /2fi-E 



q dq. 
(19) 

Resonances will appear as poles of the T matrix, 
namely as zeros of the inverse propagator of the mixed 
state. Therefore the resonance parameters are solutions 
of the equation 



A a , a (E)\ = (E-M a )S aa + G aa (E) 



= 



(20) 



with E the pole position. This equation is solved by the 
Broyden method (42l |. 

From the solution of (f2"U|) we obtain the energy and the 
total width of the resonance. However, we are faced with 
the problem of the definition of the partial width. A sim- 
ilar problem arise when one try to define the mass and 
the width of an unstable particle in a gauge independent 
way [H[ . Let assume the case of a cc bound state coupled 
to two meson states. If the two mesons are below thresh- 
old we get a mass shift of the particle mass but if they 
are above threshold the mass is also renormalized by the 
coupling and now the particle becomes unstable and ac- 



quires a width. The conventional definition of mass and 
width are in this case. 



M = M -R(0(M)), 
o 3(g(AQ) 
1 + R(0(M)) 



r 



(21) 



where Mq is the bare mass and 0(E) is the two meson 
loop. The partial width is defined by decomposing the 
numerator of the width in Eq. (|21l) into a sum of con- 
tributions of different two meson channels. However this 
cannot be done in our case because we use the complex 
value of the energy at the pole position to define the mass 
and width of the state 



E = M -0(E). 



(22) 



With the usual parametrization E = M r — iY r /2 we ob- 
tain 



r r = 2%(o(E)) 



(23) 



and we cannot follow the usual procedure to define the 
partial widths. Instead, following Ref. [HI, we start from 
the iS-matrix for an arbitrary number of cc states 



SP ,f) (E) = sf/(E) - l 2nS i (Pf - Pi) £ ^' a '(k;E)A a , a (Ey 1 r (k;E) 



(24) 



r 



where k is the on-shell momentum of the two meson state with 
and the propagator is 

A a ' a (E) = {(E- M a )S a ' a + a ' a (E)} (25) 

Expanding around the pole, which now is defined as 
\A(E)\=0: 



0' a ' a (E) = hm 5 aa (E)-0_ aa (E) 

V ; e^e E-E y ' 



A a ' a (E) - A a ' a (E) = (E-E) \5 a ' a + 0' a ' a (E) 
= (E - E)Z a ' a (E) 



Then the S'-matrix can be approximated in the neighbor- 
(26) hood of the pole as 



S^(E) = Sl'f(E) - l 2nS\P f - Pl ) £ tf* (h E)^- 



E-E 



(28) 



assuming that 



the S-matrix can be finally written as 



z a , a (E) = Y. z T^ 



(29) 
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S^^(E) = S^(E)-i2nS\P f -P i ) [^' a '(k;E)Z alx (E)- 1 / 2 



a,ot f ,A 



E-E 



(30) 



The vertex we are interested in is 

S(X C -> ff<* = ^Zx^E)- 1 ' 2 (31) 

A 

and the partial width can be defined as 

t f = [ d<Z> f \S(X c -+ f)\ 2 (32) 



where the integral is over the phase space of the final 
state with (J^nP^f = M r- 

In the case of a two meson decay Tp can be written as 

E\E 2 , 

«0/3 

(33) 



r> =2tt- 



Mr 



J2 ^ a 'Ck)Z* /x (E)- x / 2 Z(E)-^ 2 r (k) 

a', a, A 

where k^p is the onshell momentum of the two meson 
state. 

Eq. (f3"3")l does not guarantee that the sum of the partial 
widths must be equal to the total width. In fact it is 
expected that 7^ IV To solve this problem we 

define the branching ratios by [43[ 



B 



f 



and the partial widths by 

r> - BfT r . 



(34) 



(35) 



III. CALCULATIONS, RESULTS AND 
DISCUSSIONS 

Using the formalism developed in section [TT] we have 
performed a coupled channel calculation of molecular 
states and qq pairs in different sectors of the charmo- 
nium spectrum. We use the parametrization of the qq 
interaction, together with the values for quark masses 



44] 



and the strength of the 3 Pq model of Ref. 

We have analyzed the positive parity sectors from 3.8 
to 4.0 GeV where more of the new charmonium-like states 
appear. We have also calculate the effects of molecular 
structures for the J PC = 1~~ states in the region around 
4.0 GeV. 



A. Positive parity sector 

In Table U we summarized the XYZ candidates below 
4.0 GeV that we will discuss in this work together with 



State M (MeV) 


T (MeV) 


J PU 


Decay mode 


V(3872) 3871.4 ±0.6 


< 2.3 


1++ 


DD* 


X(3915) 3914±3±2 


23 ±9 


? 


D*D* 


Z(3930) 3929 ± 5 


29 ± 10 


2++ 


DD 


X(3940) 3942 ± 9 


37 ± 17 


1++ 


DD* 


V(3940) 3943 ± 17 


87 ±34 


? C = + 


J/ipu> 



Table I: Summary of candidates XYZ mesons discussed in 
this work. 



State M (MeV) J 
Xco(2P) 3909 
h c (2P) 3955 
Xci(2P) 3947 
Xc2(2P) 3968 



"0^ 

1+- 
1++ 
2++ 



Table II: Prediction from our cc model in the region around 
3970 MeV. 



the most likely J assignments and its dominant decay 
modes. In Table [TTI we show the prediction in this region 
from the model of Ref. [44[ ■ 

Although the model predicts four states in this region, 
the h c (2P) has negative C-parity and does not match 
with the data we are looking for. Then we have only 
three states predicted by the quark model whereas ex- 
perimentally one found 5 states. 

As far as masses are concern, there is one clear identifi- 
cation: the Xc2 match the -£(3930) mass (see Ref. [11] for 
more properties). For the rest of states one can briefly 
comment that the \ci is too high in mass to be V(3872). 
On the other hand the V(3940) cannot be a J p = 0++ 
state because, while a clear signal for A" (3940) — > DD* 
is seen, there is no evidence for the A (3940) in either 
DD or cuJ/ip decay channels. Then the most likely can- 
didate for A(3940) is J PC = 1++ with M = 3947 MeV. 
Concerning the A(3915), although it has not been seen 
in the DD channel due basically to its small branching 
ratio for this channel, it is a good candidate for our Xco 
state. 

With this preliminary assignments our cc model pre- 
dicts no candidates for the A(3872) and the V(3940) if 
it finally exists. 

The existence of two J PC = 1 ++ almost degenerated in 
mass, namely the A(3940) and the A(3872) suggest that 
the 1 ++ cc sector at these energies is more complicated 
than a simple cc structure. Moreover, its extremely low 
binding energy makes the A (3872) an ideal candidate to 
a DD* molecule. 
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D a D ,u 


D ± D* + 


J/ifjp J/lpLJ 




r 7 j/^ 




R2 


0.231 


-0.60 


12.40 


39.24 


7.46 


0.49 0.40 


-0.60 


8.15 


9,80 


1.20 


0.226 


-0.25 


8.00 


86.61 


4.58 


0.53 0.29 


-0.25 


5.25 


6.31 


1.20 



Table III: Binding energy (in MeV) and channel probabilities Table IV: Decay widths (in keV) of the the X(3872) states 
(in %) for the X(3872) states for two different values of the and its ratio for two different values of the 7 parameter in the 
7 parameter in the 3 Po model. 3 Po model. 



In a earlier publication |25| we performed a coupled 
channel calculation including the cc(2 3 Pi) pair together 
with the neutral and charged DD* channels. In this work 
we found a bound state with an important molecular 
component. However one can wonder if the same effects 
can be found coupling the pJ/tp and ujJ/ip neglected in 
this previous calculation. 

To elucidate this point we have performed a full cal- 
culation including the two quark cc(2 3 Pi) together with 
the D°D*°. D ± D*T : pj/i/j and LjJ/ip channels. The 
coupling of the DD* with the pJ/ip and ujJ/ip channels 
is not enough to bind the molecule and it is mandatory 
to couple the cc pair to reach a molecular bound state. 

One striking feature of the A (3872) decays is the value 
of the ratio between the A (3872) — » pJ/ip and the 
A(3872) — > utJ/tfj decay channels. As stated above this 
ratio suggests that some isospin mixing is needed to re- 
produce the experimental data. To introduce the isospin 
breaking in our calculation we will work in the charge 
basis instead of in the isospin symmetric one, allowing 
the dynamics of the system to choose the weight of the 
different components. Isospin is explicitly broken by the 
experimental meson masses. 

To have an idea of the sensitivity of the A(3872) struc- 
ture with the binding energy and having in mind that this 
ranges from 0.6 MeV to 0.25 MeV, we have fine-tuned the 
3 Po gamma parameter to get exactly this two energies. 
We can see the results in Table IIII1 

From this table one can see that the A (3872) is pre- 
dominately a D°D*° molecule with a small admixture 
(less than 10%) of the cc(2 3 Pi) state and the charged 
jj±jj*T component. The two channels pJ/4> and cjJ/tp, 
although important for the decays, are not significant 
with respect to the binding energy. 

The influence of the different components can be de- 
termine from the A(3872) decays. The ratio 



Eh 



J/ii 



-Tv°j/i> 



Ri 



R-2 



r(A(3872) 7*(25)) 
T(X(3872) ->• 7J/VO 



(36) 



is sensitive to the cc(2 3 Pi) component because in the 
molecular picture the radiative decay r(A(3872) — > 
7^(25')) is suppressed [46]. The first evidence of this 
decay was reported by BaBar with a branching frac- 

3.4 ± 1.4 [13] which suggested a 



r(X(3872)^ 7 ^(25)) _ 
r(Jf(3872)-> 7 J/V) 



tion 

rather large value of the cc(2 3 Pi) component. However, 
in 2010, using a larger sample of the B — » A(3872)A 
decay, the Belle Collaboration has not found evidences 
for the radiative decay r(A(3872) 7* (25)) giving the 



-0.60 
-0.25 



27.61 
24.18 



14.40 
10.64 



0.52 
0.44 



Table V: Strong decay widths (in keV) of the the X(3872) 
states and its ratio for two different values of the 7 parameter 
in the 3 Po model. 



Mass 


cc(2 3 A) 


D U D* [> 


D ± D* + 


J/tpp J/lpOJ 


3871.5 


8.00 


86.61 


4.58 


0.53 0.29 


3941.8 


61.09 


18.53 


16.85 


0.01 3.52 



Table VI: Masses (in MeV) and channel probabilities (in %) 
for the X(3872) and X(3940) states. 



upper limit B.TZ. < 2.1 at 90% C.L. (T^. In Table |IV] we 
show the width and the ratio for this two decays using the 
standard expression for the electric dipole transition [l]. 
One can see that the value of our cc component is small 
enough to accommodate at the experimental results. 

In Table [V] we show the calculated width for the de- 
cays X(3872) 7r+7r~J/V> andX(3872) ->■ tt+ tt~ tt° J / ^ . 
The result for the ratio is not far from the experimental 
value Pi = 0.8 ± 0.3 [9]. Although the absolute value of 
both decay widths varies with the X(3872) binding en- 
ergy their ratio is less sensitive. It is worth to notice that 
to obtain this ratio (close to the experimental value) it is 
enough to have only less than 30% of 7 = 1 component 
because, as explained in the introduction, the different 
phase spaces of the pJ/ip and the uiJ/ip channels con- 
spires with the DD* charged components to reproduce 
the experimental result. 

As stated above the formalism developed in Section |TT] 
allows us to treat simultaneously bound states and res- 
onances above and below the thresholds. When we look 
for resonances above the DD* threshold, we found a res- 
onance at M = 3941.8 MeV and T = 89.9 MeV, which 
can be identified with the A (3940). The different com- 
ponents for this resonance together with the X(3872) for 
the value 7 = 0.226 are shown in Table |VH 

As seen from Table [VI] the X(3940) has almost the 
same value of the neutral and charged DD* components. 
Therefore the isospin breaking has almost completely dis- 
appear and the resonance is basically 1 = 0. This jus- 
tifies that the coupling with the pJ/ip channel is almost 
negligible. 

This resonance decay basically through the DD* com- 
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state Mass 


r 


cc(2 J P ) 




J/tpUJ 


D S D S J/ip<j> 


X(3915) 3896.5 
Y(3940) 3970 


4.10 
189.3 


34.22 
57.27 


46.67 
35.32 


9.42 
0.15 


9.67 0.03 
5.72 1.54 



Table VII: Mass and total width (in MeV) and channel prob- 
abilities (in %) for the X(3915) and Y(3949). 



ponent being the branching ratio for the different chan- 
nels S.ft.(X(3872) -> DD*) = 0.89, £."K.(X(3872) -> 
uJ/^i) = 0.1 and £."K.(A(3872) -> pJ/V») = 3 • 10~ 4 . 
In this way the puzzle between this two states with the 
same quantum numbers seems to be solved in a satisfac- 
tory way. 

The situation in the J PC = ++ is also puzzling. Be- 
fore the discovery of the A(3915) signal, it was supposed 
that the F(3940) could be the J PC = 1++ or J PC = 0++ 
state. If we look to the Table HH our cc model predicts 
a J PC = 0++ state at 3909 MeV and so the F(3940) 
seems to be too high. Despite the fact that one of the 
states is narrow (the A" (3915)) and the other is broad (the 
Y(3940)), lately there has been a tendency to consider 
that both are the same state. We have performed the 
same calculation as before but for the J PC — ++ sec- 
tor. We include the cc(2 3 Pq) channel together with the 
molecular channels DD (3736.05 MeV), uJ/i[> (3879.56 
MeV), D S D S (3936.97 MeV), 0J/V> (4116.0 MeV) where 
the channel thresholds are given between parenthesis. 
The results are shown in Table IVTTl 

We find a narrow state which can be identified with 
the X(3915). This resonance is basically a mixture of 
cc and a DD molecular component. Moreover a second 
wide resonance appears at M = 3970 MeV. This reso- 
nance can be the old F(3940), today disappeared from 
the Particle Data Group. Its dominant component is the 
cc(2 3 Po) although the contribution of the molecular DD 
one is also significant. This structure provide a possible 
explanation for the unusual decay mode uJ/tp through 
the rescattering ++ — > DD —> uJ/ip. 

Then our calculation shows that the X(3915) and the 
y(3940) may be two different resonances as measured by 
Belle. 



B. Negative parity sector 

In the J PC = sector one can find a similar sit- 
uation to those studied in the previous sector. In fact 
the '0(4040) resonance with a mass just above the D*D* 
threshold has been proposed long ago [11] as a candidate 
to a molecular state. To asses this possibility we have 
performed a coupled channel calculation including the 
cc(3 3 S'i) and cc(2 3 L>i) states with masses M = 4097.615 
MeV and M = 4152.715 MeV respectively, together with 
the channels DD, DD* , D*D* , D S D S , D S D* S and D*D *. 
The results of the calculation are shown in Table IVIIIl 

The first striking outcome of the calculation (see Ta- 



ble IVIIIf) is the appearance of a narrow state at M = 
3994.6 MeV. A state with this characteristics has been 
reported by BaBar [4!| in the study of exclusive initial- 
state-radiation production of the DD system. Its exper- 
imental mass and width are respectively M — 3943 ± 
17 ± 12 MeV and T = 52 ± 8 ± 7 MeV. This resonance 
has also been seen by the Belle Collaboration [5(|. The 
second significant result is that, due to the molecular 
mixing, the cc(2 3 £>i) state becomes the most important 
component of the ■0(4040) and not the cc(3 3 Si) as usu- 
ally attributed by the naive quark model. Moreover the 
0(4040) acquires a significant DD* (23.49%) and D*D* 
(25.81%) components while the 0(4160) is predominantly 
a cc(3 3 5i) state with small contributions (< 20%) of dif- 
ferent DD and D S D S channels. 

These new assignments have an important influence 
on the decay ratios of these two resonances. It is well 
known that the ratios of branching fractions involving 
these two resonances shows significant discrepancies with 
model predictions, specially with the 3 Po model. The 
predictions of our calculation are shown in Table IIXI to- 
gether with two different quark models. As seen, the 
new assignment improve the overall agreement with the 
experimental data although some discrepancies remains. 

Another consequence of the important molecular com- 
ponent of the '0(4040) is that it can give rise to an 
enhancement of some specific decay channels like the 
V>(4040) -> 77J/0 through the presence of {D*D*) s =o 
pairs in the ^(4040) internal structure. This mechanism 
has been recently proposed by Voloshin [H3| as a tool to 
identified possible molecular structures in this states. 



IV. SUMMARY 

During the years charmonium spectrum below DD 
threshold was a well described system in quark mod- 
els. However since the last 10 years new states above 
this threshold has been measured with properties which 
can be hardly described in naive quark models. Obvi- 
ously at these energies threshold effects has to be taken 
into account which sometimes has been referred as an 
unquenching of the quark model. 

In this paper we address this issue focusing on the pos- 
sible influence of the molecular structures on the charmo- 
nium spectrum. In the framework of a constituent quark 
model, we have developed a formalism which allows us to 
coupled cc states with four quark molecular states below 
and above the different thresholds. We study the positive 
parity sector in the mass region of the X(3872) and the 
negative parity sector around masses of 4.0 GeV. 

We describe the A (3872) resonance as a 

JPC = 1+ + 

mixture of neutral and charged DD* molecular states 
and a less than 10% contribution of the cc(2 3 Pi) state 
what however is enough to describe the electromagnetic 
decays of the resonance. The isospin breaking showed by 
the data is also well explained with this configuration. 
Together with this resonance in the 1 ++ sector appears a 
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M{MeV) 


cc3 3 5j 


cc2 a Di 


DD 


DD* 


D*D* 


d s d s d s d; d* 3 d* 


3994.6 -ill.60 


31.56 


3.0 


2.49 


36.44 


117.75 


7.53 0.523 0.71 


4048.4 - i7.54 


0.92 


36.15 


2.99 


23.49 


25.81 


8.86 0.924 0.85 


4123.9 - i71.ll 


59.01 


0.98 


2.13 


6.84 


19.119 


0.75 3.37 7.73 



Table VIII: Mass (in MeV) and channel probabilities (in %) for the J = 1 sector in the 4.0 GeV region. 





ratio 




Measurements 


J Po f51] 


C A [52] This work 


B.7^(0(4O4O) 


-> DD)/B.K. (0(4040) - 


> DD*) 


0.24 ±0.05 ±0.12 


0.003 


0.0003 


0.07 


5.^.(^(4040) - 


V D*D*)/B.TZ.{ip(4040) 


-> DD*) 


0.18 ±0.14 ±0.03 


1.0 


1.0 


0.61 


S.ft.(V>(4160) - 


-> DD)/B.H.(tI>(4160) - 


• D*D*) 


0.02 ±0.03 ±0.02 


0.46 


0.008 


0.0.5 


Z3.ft.(V>(4160) - 


> DD*)/B.K.(ip(<il60) - 


» D*D*) 


0.34 ±0.14 ±0.05 


0.011 


0.16 


0.08 



Table IX: Ratios of branching fractions for the two ip resonances. Experimental data are from ref [531 ]. 



second state whose properties are compatible with those 
of the A (3940) state. 

We found two resonances with J PC = ++ quantum 
numbers. The first one, with important cc and DD 
components, may be identified with the A(3915). We 
found also a second broad resonance in the mass region 
M = 3940 MeV which could be assign to the signal 
seen by the Belle Collaboration and which for sometime 
was called F(3940), although recently, probably due to 
the lack of quantum numbers to accommodate this reso- 
nance, was subsummed under the name A(3915). 

Concerning the negative parity sector, we confirm the 
old suggestion of De Rujula et al. [H| that the 0(4040) 
resonance is mostly a molecular state. Moreover the 
coupling with these molecular structures change the cc. 
quantum numbers of the 0(4160) which acquires an im- 
portant 3 Si component, contrary to the usual hypothesis 



that is a 3 I?i state. This new assignment has important 
consequences on the decay branching ratios. Finally the 
molecular components of these two resonances open the 
possibility of enhance the probability of new decay chan- 
nels to detectable levels which deserves further studies. 
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